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1. Introduction 

Electric dipole spin resonance (EDSR) is a method to electrically manipulate electron 
spins. In this technique, two fields are applied; a static magnetic field B and an 
oscillating electric field E(t) resonant with the electron precession (Larmor) frequency [U 
El [3l Hj. Spin resonance techniques are of interest for quantum computing schemes 
based on single electron spins, because they allow arbitrary one-qubit operations [5]. 
Single-spin EDSR is a particularly desirable experimental tool because it allows spin 
manipulation without time-dependent magnetic fields, which are difficult to generate 
and localize at the nanoscale j6j EJ E] . 

Achieving EDSR requires a mechanism to couple E to the electron spin <r. This 
coupling can be achieved by the traditional spin-orbit interaction, which couples a to 
the electron momentum k, or by an inhomogeneous Zeeman interaction, which couples 
a to the electron coordinate r [10l HI HU 021 E]- Single-spin EDSR has recently been 
achieved in quantum dots using both techniques [IU [15] . 

Recently, we presented an experimental and theoretical study of a novel EDSR effect 
mediated by the spatial inhomogeneity of the hyperfine nuclear field [16]. An electron 
moving under the influence of the electric field E(t) experiences this inhomogeneity 
as an oscillating hyperfine coupling which drives spin transitions. In this paper, 
we illuminate the underlying physics and present new experimental data on a still 
unexplained phenomenon at half the resonant frequency. 

This EDSR effect is observed via spin-blocked transitions in a few-electron GaAs 
double quantum dot [T7]. As expected for a hyperfine mechanism, but in contrast to the 
k — cr-coupling mediated EDSR, the resonance strength is independent of B at low field 
and shows, when averaged over nuclear configurations, no Rabi oscillations as a function 
of time. We find that at large B driving the resonance creates a nuclear polarization, 
which we interpret as the backaction of EDSR on the nuclei [THl EE EE EOj [21]. Finally, 
we demonstrate that spins can be individually addressed in each dot by creating a local 
field gradient. 

2. Device and measurement 

The device for which most data is presented (Figure 1(a)) was fabricated on a 
GaAs/Al .3Ga .7As heterostructure with two-dimensional electron gas (2DEG) of 
density 2 x 10 15 m -2 and mobility 20 m 2 /Vs located 110 nm below the surface. Voltages 
applied to Ti/Au top gates locally deplete the 2DEG, defining a few-electron double 
quantum dot. A nearby charge sensing quantum point contact (QPC) is sensitive to 
the electron occupation (N Ll N R ) of the left (N L ) and right (N R ) dots [221 123]. The 
voltages Vl and Vr on gates L and R can be rapidly pulsed; in addition, L is coupled 
to a microwave source. The static magnetic field B was applied in the plane of the 
heterostructure, and measurements were performed in a dilution refrigerator at 150 mK 
electron temperature. 
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Figure 1. (a) Micrograph of a device lithographically identical to the one measured, 
with schematic of the measurement circuit. Quantum dot locations are shown by 
dashed circles, and a bias V s a drives sequential tunneling in the direction marked 
by black arrows. The conductance g s of the QPC on the right is sensitive to the 
dot occupation. The direction of the magnetic field B and the crystal axes are 
indicated, (b) QPC conductance g s measured at V s d ~ 600 fieV near the (l,l)-(0,2) 
transition. Equilibrium occupations for different gate voltages are shown, as are gate 
voltage configurations during the measurement/reinitialization (M) and manipulation 
(C) pulses. The two white dashed triangles outline regions where transport is not 
Coulomb blocked; the solid black line outlines where spin blockade is active. A plane 
background has been subtracted, (c) Energy levels of the double dot during the pulse 
cycle (See text). 



The characteristic feature of tunnel-coupled quantum dots is a discrete electron 
energy spectrum. An overall shift to the spectrum, proportional to the electron 
occupation, is induced by Vl and Vr, which therefore determine which occupation is 
energetically favoured. Figure 1(b) shows the QPC conductance q§ as a function of 
Vl and Vr; different conductances correspond to different (N L ,N R ). For most Vl,Vr 
configurations, only one value of (N L ,N R ) is energetically accessible; these correspond 
in Figure 1(b) to regions of uniform g s . 

A bias V s d applied across the device drives electron transport via sequential 
tunneling subject to two constraints [24J. The first constraint, Coulomb blockade, arises 
because for most gate configurations electrostatic repulsion prevents additional electrons 
from tunneling onto either dot. This constraint inhibits transport except when V^, Vr are 
tuned so that three occupation configurations are near-degenerate. The energy cost of an 
extra electron tunneling through the device is then small enough to be provided by the 
bias voltage. Values of Vl and Vr satisfying this condition correspond to the two white 
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dashed triangular regions marked in Figure 1(b), for which transport is permitted via the 
transition sequences (0, 2) (0, 1) (1, 1) (0, 2) or (0, 2) (1, 2) (1, 1) (0, 2). 

A second constraint, spin blockade, is caused by the Pauli exclusion principle, which 
leads to an intra-dot exchange energy J 02 in the right dot [25, 26j. As shown in the first 
panel of Figure 1(c), the effect of this exchange is to make the (1, 1) —> (0, 2) transition 
selective in the two-electron spin state, inhibited for triplet states but allowed for the 
singlet. Although the hyperfine field difference between dots rapidly converts the m s = 
component T of the blocked triplet T to an unblocked singlet S, decay of m s = ±1 
components T± requires a spin flip and therefore proceeds much more slowly. This spin 
flip becomes the rate-limiting step in transport, and so the time- averaged occupation 
is dominated by the (1,1) portion of the transport sequence [26J. Gate configurations 
where spin blockade applies correspond to the black solid outlined region of Figure 1(b); 
inside this region, g s has the value corresponding to (1,1). Any process that induces 
spin flips will partially break spin blockade and lead to a decrease in g s . 

Unless stated otherwise, EDSR is detected via changes in g s while the following 
cycle of voltage pulses Vl and Vr [9] is applied to L and R (Figure 1(c)). The cycle 
begins inside the spin blockade region (M in Figure 1(b)), so that the two-electron state 
is initialized to (1, 1)X± with high probability. A ^1 /xs pulse to point C prevents electron 
tunneling regardless of spin state. Towards the end of this pulse, a microwave burst of 
duration tedsr at frequency / is applied to gate L. Finally the system is brought back 
to M for ^3 /is for readout/reinitialization. If and only if a spin (on either dot) was 
flipped during the pulse, the transition (1,1) — ► (0,2) occurs, leading to a change in 
average occupation and in g s . If this transition occurs, subsequent electron transitions 
reinitialize the state to (1, 1)T± by the end of this step, after which the pulse cycle is 
repeated. This pulsed EDSR scheme has the advantage of separating spin manipulation 
from readout. 

Changes in g s are monitored via the voltage Vqpc across the QPC sensor biased at 
5 nA. For increased sensitivity, the microwaves are chopped at 227 Hz and the change in 
voltage S Vqpc is synchronously detected using a lock-in amplifier. We interpret 5 Vqpc 
as proportional to the spin-flip probability during a microwave burst, averaged over the 
100 ms lock-in time constant. 

3. EDSR spectroscopy 

Resonant response is seen clearly as B and / are varied for constant tedsr = 1 //s 
(Figure 2.) A peak in 5 Vqpc, corresponding to a spin transition, is seen at a frequency 
proportional to B. This is the key signature of spin resonance. From the slope of the 
resonant line in Figure 2 a ^-factor \g\ = 0.39 ± O.Olis found, typical of similar GaAs 
devices [271 [28]. We attribute fluctuations of the resonance frequency (inset of Figure 2) 
to Overhauser shift caused by the time-varying hyperfine field acting on the electron 
spin. Their range is ~ ±22 MHz, corresponding to a field of ~ 4 mT, consistent with 
Overhauser fields in similar devices [29. 1301 13T] . 
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Figure 2. Signal of spin resonance SVqpc as a function of magnetic field B and 
microwave frequency /. EDSR induces a breaking of spin blockade, which appears 
as a peak in the voltage across the charge sensor SVqpc at the Larmor frequency. 
Field- and frequency-independent backgrounds have been subtracted. Inset: Jitter of 
resonant frequency due to random Overhauser shifts. 



Information about the EDSR mechanism can be obtained by studying the peak 
height as a function of duration, strength, and frequency of the microwave burst 
(Figure 3). To reduce the effects of the shifting Overhauser field, the microwave source 
is frequency modulated at 3 kHz in a sawtooth pattern with depth 36 MHz about a 
central frequency /. The resonance line as a function of tedsr is shown in the inset of 
Figure 3(a). For equal microwave power at two different frequencies /, the peak heights 
SVqPq are plotted in Figure 3(a) (main panel). The two data sets are similar in turn-on 
time and saturation value; this is the case for frequencies up to / = 6 GHz. From similar 
data (insets of Figure 3(b)), using theory to be described, we extract the dependence of 
the spin-flip rate Qr on microwave power Pmw shown in the main panel of Figure 3(b). 
Coherent Rabi-type oscillations in 5Vqp^(tedsr) are not observed for any microwave 
power or magnetic field over the range measured. 

The ^-independence of the EDSR strength rules out spin-orbit mediated EDSR 
of the k — a type (either Dresselhaus or Rashba), for which the Rabi frequency is 
proportional to B [H CH [H] . This is in contrast to the results of [14] , where the spin- 
orbit effect was found to dominate in a similar device to ours. A possible explanation 
is the device orientation relative to B and the crystal axes. In both our experiment 
and [14J, the gate geometry suggests a dominant E(t) oriented along one of the diagonal 
axes ([110] or [110]), leading to an in-plane spin-orbit effective field B^p perpendicular 
to E(t). In our geometry (see Figure 1(a)), this orientation of B^j? is parallel to B, 
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Figure 3. (a) Measured EDSR peak strength SVq^ (symbols) versus microwave 
pulse duration tedsr f° r two frequencies at equal power, along with theoretical fits 
(curves) obtained by numerically evaluating and scaling Equation (4) (see text). Inset: 
Raw data from which the points in the main figure are extracted. Each vertical cut 
corresponds to one point in the main figure. Jitter in the field position of the resonance 
reflects time-dependent Overhauser shifts, (b) Spin-flip rate CIr as a function of applied 
microwave power Pmw, along with a fit to the form oc VPmw (dashed line). Insets: 
^Vqpc versus tedsr fc> r two values of the microwave power, showing the fits from which 
points in the main figure are derived. 



and therefore ineffective at driving spin transitions. In the geometry of [14], B is 
perpendicular to B^j?, so that the k — a spin-orbit mechanism becomes more efficient . 

4. Theory 

A theoretical description of 5Vqp^(Yedsr) and its dependence on B and Pmw can be 
obtained by modeling EDSR as arising from the coupling of an electron in a single dot to 
an oscillating electric field E(t) and the hyperfine field of an ensemble of nuclei [|| [T61I32] . 
Then the center of the dot oscillates as R(i) = — e~E(t) /itiuJq, where m is the electron 
effective mass, and uj is its confinement frequency in a parabolic dot. As a result, the 

£ There exists some physical similarity between the hyperfine mechanism of EDSR described in this 
paper and EDSR due to the coupling of electron spin to a random exchange field in semimagnetic 
semiconductors [33] . 
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Hamiltonian of the hyperfine coupling of the electron spin S = a/2 to nuclear spins Ij 
becomes time dependent, = A£j5(r + R(t) — Tj)(lj • S). Here A is the hyperfine 
coupling constant and the summation over j runs over all nuclear spins. After expanding 
i/hf in R(t) (assumed small compared to the dot size) and averaging over the orbital 
ground-state wave function ^o( r ) of the dot, the time dependent part of becomes 
Hhf(t) = J(t) • <r, where J(t) is an operator in all Ij. Choosing the z-axis in spin space 
along B, the components of J(t) are J z = \A J^ . ^(r^ij and 

J ±W = ^2 E^(ri)E(t) • VVbfo)^. (1) 



J 

The time-dependent off-diagonal components J±(t) drive EDSR, while the quasi- 
static diagonal component J z describes detuning of EDSR from the Larmor frequency 
ujl by an amount uj z randomly distributed as p(uj z ) = exp(— cj|/A 2 )/(Ay / 7r) [34J. 
The dispersions A of the detuning and Qr of the Rabi frequency are the root-mean- 
square values of J z and J± respectively. Whereas J z is dominated by fluctuations of Ij 
symmetric about the dot centre, J± is dominated by fluctuations antisymmetric in the 
E direction because E • V^o(r) is odd with respect to the E projection of r. Finally, 



A I(I+l)mcu no n eEA I(I + l)n 
A - 2ft V 2M ' Qr ~ JKJ-J 8nd ' (2) 
with / = 3/2 for GaAs, no the nuclear concentration, and d the vertical confinement. 
It is seen that Qr is independent of B\ this is in contrast to EDSR mediated by the 
conventional k — a spin-orbit coupling, where Kramers' theorem requires that the Rabi 
frequency vanish linearly as B — > [H [35j [[3] . 

In an instantaneous nuclear spin configuration with detuning 8uj = 2nf — (ujl +uj z ) 
and Rabi frequency the spin-flip probability from an initial | spin state is [36] : 

a? 



^ (TEDSR) = (^fT^ sin2 



(5u/2f + TEDSR 



(3) 



(We neglect the electron spin relaxation and nuclear-spin dynamics, which remain slow 
compared with the Rabi frequency even in the EDSR regime [3T1 [32].) To compare 
with the time-averaged data of Figure 3, we average Equation ^ over cu z with weight 
p((jj z ) and over Vt with weight p{£l) = 2Qexp(— Vt 2 / Vt 2 ^ / Vt 2 R . This latter distribution 
arises because the J± acquire Gaussian-distributed contributions from both J? and I y - 
components of the nuclear spins, hence it is two-dimensional. Averaging over uu z and Q 
results in a mean-field theory of the hyperfine-mediated EDSR. The resulting spin-flip 
probability 

/+oo r+oc 
du z p(u z ) / dCl p(Cl)p^ (tedsr) (4) 
oo JO 

shows only a remnant of Rabi oscillations as a weak overshoot at tedsr ~ ^Ir 1 - 
The absence of Rabi oscillations is characteristic of hyperfine-driven EDSR when 
the measurement integration time exceeds the nuclear evolution time [37J, and arises 
because J± average to zero. 
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4-1- Comparison with data 

To compare theory and experiment, the probability ^(tedsr! A, Qr) is scaled by a QPC 
sensitivity Vq PC to convert to a voltage SVq^. After scaling, numerical evaluation of 
Equation' (4) gives the theoretical curves shown in Figure 3(a). The parameters that 
determine these curves are as follows: The Larmor frequency spread, A = 2n x 28 MHz, 
is taken as the quadrature sum of the jitter amplitude seen in Figure 2 and half the 
frequency modulation depth, whereas FIr and Vq PC are numerical fit parameters. The 
44 mT data (green curve in Figure 3(a)) give Qr = 1.7 x 10 6 s _1 and Vq PC = 2.4 /xV. 
Holding Vq PC to this value, the 550 mT data give Qr = 1.8 x 10 6 s _1 (blue curve in 
Figure 3(a)) and the 185 mT data give the dependence of Qr on microwave power Pmw 
shown in Figure 3(b). The Rabi frequency Qr increases as \APmw (Figure 3(b)) and is 
independent of £>, both consistent with Equation (1). The ^-independence of Qr — also 
evident in the EDSR intensity in Figure 2 — and the absence of Rabi oscillations support 
our interpretation of hyperfine-mediated EDSR in the parameter range investigated. 

Estimating Hujq ~ 1 meV [28], E ~ 6 x 10 3 Vm" 1 at maximum applied power [§J 
d ~ 5 nm, and using values from the literature n = 4 x 10 28 m -3 and An =90 /ieV 
[38] we calculate Qr ^ 11 x 10 6 s _1 , an order of magnitude larger than measured. The 
discrepancy may reflect uncertainty in our estimate of E. 

Above, we generalized a mean-field description of the hyperfine interaction [39. 13^] 
to the resonance regime. Justification for this procedure was provided recently in [32J. 
A distinctive feature of the mean-field theory is a weak overshoot, about 10 - 15%, that 
is expected in the data of Fig. 3(a) before 5Vq P ^(tedsr) reaches its asymptotic value at 
7" edsr —> °°- No overshoot is observed in the 550 mT data (blue symbols in Figure 3(a)), 
which was taken in a parameter range where an instability of the nuclear polarization 
begins to develop; see Section 5. For the 44 mT data (green symbols in Figure 3(a)), a 
considerable spread of experimental points does not allow a specific conclusion regarding 
the presence or absence of an overshoot. The theory of [32] suggests that the existence 
of the overshoot is a quite general property of the mean-field theory. However, after 
passing the maximum, the signal decays to its saturation value vary fast, with Gaussian 
exponent e r t edsr m By contrast, the first correction to the mean-field theory decays 
slowly, as 1/ (^^r t edsr) > where is the number of nuclei in the dot. As a result, the two 
terms become comparable at tedsr ^ \/ln N /Or, which should make the maximum less 
pronounced. Because for N ~ 10 5 the factor \/hiN ~ 3, the corrections to the mean- 
field theory manifest themselves surprisingly early, at times only about tedsr ~ 3/Ar, 
making the overshoot difficult to observe. 
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Figure 4. (a) Shift of the resonance frequency with time at constant B — 2390 mT, 
showing build-up of nuclear polarization over ~ 200 s. (b) A scheme to allow larger 
polarizations: the microwave frequency is repeatedly scanned over the resonance 
while B is swept upwards. Nuclear polarization partly counteracts £?, moving the 
resonance away from its equilibrium position (black diagonal line) by up to 840 mT. (c) 
Similar data taken at lower microwave power and opposite frequency sweep direction, 
showing approximately the equilibrium resonance position. (Grey scale as in (b)). (d) 
Similar data as in (b), with faster sweep rate, showing more clearly the displacement 
and subsequent return to equilibrium of the resonance. marks the escape of the 
resonance from the swept frequency window. In all plots, arrows denote frequency 
sweep direction. 



5. Nuclear polarization 

Consistent with a hyperfine mechanism, this EDSR effect can create a non-equilibrium 
nuclear polarization [20] . If / is scanned repeatedly over the resonance at high power, a 
shift of the resonance develops (Figure 4(a)), corresponding to a nuclear spin alignment 
parallel to B. The effect is stronger at higher £>, and saturates over a timescale ~ 200 s. 
In Figure 4(b), we show how to build up a substantial polarization: While slowly 
increasing £>, we scan / repeatedly downwards, i.e., in the direction which tracks the 
moving resonance. The resonance frequency remains approximately fixed, showing that 
the developing polarization compensates the increase in B. From the maximum line 
displacement from equilibrium, an effective hyperfine field of 840 mT can be read off, 
corresponding to a nuclear polarization of ~ 16%. Figure 4(c) shows similar data 

§ The maximum power is limited by non-resonant lifting of spin blockade, which we take to indicate a 
microwave amplitude exceeding the gate voltage from C to the nearest charge transition. The data in 
Figure 3(a) and the last data point in Figure 3(b) use power 2±1 dB below this threshold, corresponding 
to 3.2 mV. Dropped uniformly across the whole device this voltage gives a field E ~ 3 x 10 3 Vm" 1 . 
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Figure 5. (a) A device similar to that of Figure 1, incorporating a micromagnet. (b) 
Total field magnitude B tot (right axis) and the x and z components of the micromagnet 
contribution B mag (left axis), simulated at y = for external field B = 200 mT along 
z (out of the plane). £?™ ag vanishes by symmetry. The gate layout is shown in the 
background, (c) The associated split EDSR line. The lower resonance is stronger, as 
expected if the left electron is confined close to the minimum of B tot 



for lower power and opposite frequency sweep direction, indicating the approximate 
equilibrium line position. Figure 4(d), similar to Figure 4(b) but with a faster sweep 
rate, makes the displacement and eventual escape of the resonance clearer although the 
maximum polarization is less. 

The resonance shift is observed to be towards lower frequency, corresponding 
to a nuclear polarization parallel to B. This can be understood if the pulse cycle 
preferentially prepares the electron ground state T + over T_, either because it is more 
efficiently loaded or because of electron spin relaxation. EDSR then transfers this 
electron polarization to the nuclei [21]. We emphasize that the line shift is opposite 
to what is given by the usual Overhauser mechanism for inducing nuclear polarization 
via electron resonance [4(2 [18] . 

6. Addressing individual spins 

In quantum information applications, it is desirable to address individual spins 
selectively [5]. A scheme to allow this is presented in Figure 5. In an otherwise similar 
device (Figure 5(a)), we incorporated a 100 nm thick micron-scale permalloy (84% Ni, 
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Figure 6. Spin resonance signal (measured in conductance) in the device of Figure 
5(a). The EDSR signal shows up as a decrease in conductance as expected at frequency 
corresponding to\g\ = 0.45 (marked with dashed line.) An additional signal of opposite 
sign appears at exactly half this frequency (dotted line). As in Figure 2, field- and 
frequency-independent backgrounds have been subtracted. The horizontal features at 
0.5 and 1.5 GHz result from resonances of the microwave circuit. 



16% Fe) magnet over 35 nm of atomic-layer-deposited alumina [T2l|4T]. This device was 
measured with external field B normal to the heterostructure plane. A finite-element 
simulation of the field B mag due to the micromagnet, assuming complete permalloy 
magnetization along B, yields the field profiles shown in Figure 5(b). The difference 
in total field B tot = |B + B mag | between dots is - 7 mT. As expected, the EDSR line 
measured in this device is frequently split (Figure 5(c)). The splitting, 10 — 20 mT 
depending on precise gate voltage and pulse parameters, is not observed without the 
magnet and presumably reflects the field difference between dots. Since this splitting is 
considerably larger than the Overhauser field fluctuations, spins in left and right dots 
can be separately addressed by matching / to the local resonance condition fT5] . 

7. Open issues and discussion 

Finally, we discuss unexplained behavior observed only in the device of Figure 5(a). For 
the data described in this section, a simplified measurement scheme is used: Rather 
than applying gate pulses, the device is configured in the spin blockade region (point M 
in Figure 1(a)) throughout. Microwaves are applied continuously, and spin resonance is 
detected by directly measuring the QPC conductance g s . 

As well as the EDSR signal at full frequency / = g/LisB/h^ an unexpected half- 
frequency signal is sometimes seen (Figure 6.) Furthermore, depending on the exact gate 
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configuration, both full-frequency and half-frequency signals can have either sign; the 
change in g s at full frequency is usually negative as expected, but sometimes positive 
close to degeneracy of (1,1) and (0,2) charge configurations, where spin blockade is 
weakest [29]; by contrast, the change in g s at half frequency is usually positive but 
sometimes negative far from degeneracy. For most gate configurations, full-frequency 
and half-frequency signals have opposite sign, as seen in Figure 6. 

A half- frequency response is as far as we know unprecedented in spin resonance, 
and suggests second harmonic generation (SHG) from the microwave field. SHG is 
generally a non-linear phenomenon; it occurs for example in optical materials with non- 
linear polarizability [12] and in non-linear electronic components such as diodes. In 
our system, a hyperfine field at a harmonic of the microwave frequency arises if the 
confinement potential is non-parabolic. 

However, SHG alone does not explain the sign of the conductance change seen 
at half-frequency in Figure 6. A positive signal could in principle be caused by an 
admixture of the (0,1) charge state; but it is observed even for the gate configurations 
where (0,1) is energetically inaccessible (in the top right of the spin blockade region of 
Figure 1(b)). Also, there is no reason why (0,1) should be admixed for one resonance 
but not the other. These anomalous behaviours are therefore left unexplained. 
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